c-axis charge gap and its critical point in the heavily doped Ba(Fei_^Co2:)2As2 
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Temperature-dependent inter- plane resistivity, pciT), was used to characterise the normal state 
of the iron-arsenide superconductor Ba(Fei_2,Coa;)2 As2 over a broad doping range < a; < 0.50. 
The data were compared with in-plane resistivity, pa{T), and magnetic susceptibility, taken 
in HJ-C, as well as Co NMR Knight shift, ^^K, and spin relaxation rate, 1/T\T. The inter-plane 
resistivity data show a clear correlation with the NMR Knight shift, assigned to the formation 
of the pseudo-gap. Evolution of Pc{T) with doping reveals two characteristic energy scales. The 
temperature of the cross-over from non-metallic, increasing on cooling, behavior of pc{T) at high- 
temperatures to metallic behavior at low temperatures, T* , correlates well with an anomaly in all 
three magnetic measurements. This characteristic temperature, equal to approximately 200 K in 
the parent compound, x=Q, decreases with doping and vanishes near x* f»0.25. For doping levels 
X > 0.166, an additional feature appears above T* , with metallic behavior of Pc{T) found above the 
low-temperature resistivity increase. The characteristic temperature of this charge-gap formation, 
TcG, vanishes at xqg —0.30, paving the way to metallic, T-linear, pc(T) close to xqg and super- 
linear T-dependence for x > xcG- None of these features are evident in the in-plane resistivity 
pair). For doping levels x < xca, xiT) shows a known, anomalous, T-linear dependence, which 
disappears for x > xca- These features are consistent with the existence of a uniaxial charge gap, 
accompanying formation of the magnetic pseudogap, and its critical suppression with doping. The 
inferred c-axis charge gap reflects the three-dimensional character of the electronic structure and of 
the magnetism in the iron arsenides. 

PACS numbers: 74.70.Dd,72.15.-v,74.25.Jb 
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The metallic state of the, until recently, only known 
high temperature superconductors, the compounds based 
of Cu-0 elements and frequently referred to as the 
cuprates is characterized by a plethora of anomalies. 
At low doping levels, anomalous behaviours are found 
in the temperature-dependent resistivity, magnetization, 
NMR Knight shift and relaxation rate, as well as in spec- 
troscopic data 0. These behaviours are consistent with 
a decrease in the density of states at low temperatures, 
usually assigned with pseudogap formation. The phe- 
nomenology and fc-space distribution of the pseudogap 
in the cuprates is now well established {sj, however, its 
microscopic origin is still debated [^]. Main theories and 
experiments link it to two neighboring phases, an antifer- 
romagnetic Mott-insulator, with pseudogap arising due 
to exotic magnetism [5], or to a superconducting phase, 
as an effect of the preformed superconducting pairs @. 
The pseudogap is universally observed in both hole and 
electron Q doped cuprates, though it is much more pro- 
nounced in the former. 

Discovery of superconductivity with high critical tem- 
peratures in FeAs-based materials Q, breaking the 
monopoly of the cuprates, naturally raises the question 
about the common features of the two families jQ]. It 
fuels the hopes that one day the enigmatic mechanism of 
high temperature superconductivity will be understood. 
One of the important features to understand from such 
comparison, is a possible link between superconductivity 
and the pseudogap. 



Features consistent with pseu dogap are indeed ob- 
served in the hole doped RFeAsO |10l4l4| (R= rare earth, 
1111 compounds in the following). A clearly decreased 
density of states is found in ARPES measurements in 
(Ba,K)Fe2As2 (BaK122 compounds in the following) fisj . 
Because the parent compounds of iron pnictides are met- 
als, the pseudogap here is believed to arise from nesting 
instability |1(|. 

On the contrary, the experimental situation in electron 
doped Ba(Fei_2;Co2;)2As2 (BaCol22 in the following) is 
less clear. NMR studies suggest the existence of a pseu- 
dogap over the whole doping range, from magnetically 
ordered parent compound to overdoped superconductor, 
with a characteristic temperature of 560 K ±150 K at op- 
timal doping 17[; ARPES found a tiny feature just above 



the superconducting Tc 18|, whereas the in-plane resis- 



tivity does not reveal any pseudogap-like features [19| 
and is well described in a broad composition range by a 
sum of T-linear and T^-contrinutions [20| . 

We have recently undertaken extensive anisotropic 
electrical resistivity measurements on parent and opti- 
mally doped Ba(Fei_2:Coa;)2As2 2l|-(23|. In addition to 
a small ac-anisotropy, we found different temperature de- 
pendencies of the in-plane and inter-plane electrical re- 
sistivity. Here we report a systematic study of the evolu- 
tion of the inter-plane resistivity with doping. We show 
that the anomalies in the inter-plane resistivity reflect the 
existence of the enigmatic pseudogap state in BaCol22. 
Clear correlation with NMR measurements in BaCol22 
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as a function of doping, [2J] and the lack of any associ- 
ated features in the in-plane transport, suggest uniaxial 
symmetry of the pseudogap. 

Tracking the evolution of the characteristic features 
of the temperature-dependent inter-plane resistivity with 
doping we found a critical concentration, xqg ~0.30, be- 
yond which the pseudogap features disapper. Our mag- 
netization measurements show that this corresponds to 
the concentration at which the magnetic susceptibility 
loses its anomalous T-linear increase at high tempera- 
tures. At the critical concentration, the Pc{T) is very 
close to linear. This evolution of the inter-plane electri 
cal resistivity suggests a (quantum) critical point [IT 
on the edge of the pseudogap state. 



EXPERIMENTAL 

Single crystals of BaFe2 As2 doped with Co were grown 
from a starting load of metallic Ba, FeAs and CoAs, as 
described in detail elsewhere [13] • Crystals were thick 
platelets with sizes as big as 12x8x 1 mm"^ and large faces 
corresponding to the tetragonal (001) plane. The actual 
content of Co in the crystals was determined with wave- 
length dispersive electron probe microanalysis and is the 
x- value used throughout this text. 

In our study of resistivity anisotropy in 
Ba(Fei„a;Coa;)2As2, undoped a;=0 [22| and opti- 
mally doped x=0.074 



2l|, we have found that special 



care must be taken for measurements in configurations 
with current along the tetragonal c-axis so as to avoid 
effects associated with the exfoliation of the samples. 
Cutting and shaping into transport samples inevitably 
introduces cracks, which affect the effective geometric 
factors of the samples and, in case the cracks are 
deep, can produce admixture of the in-plane resistivity 
component. A strong tendency to exfoliate prevents the 
cutting of samples with c 3> a. This limitation puts 
severe constraints on the measurement technique. 

Samples for electrical resistivity measurements with 
current flow along the tetragonal c axis {pc) were cut 
into (0.3-0.7) X (0.3-0.7) X (0.1-0. 5)mm3 (a x 6 x c) slabs. 
All sample dimensions were measured with an accuracy 
of about 10%. Contacts to the samples were made by 
attaching silver wires using ultrapure tin, resulting in 
an ultra low contact resistance (less than 10 /Ltfi) [26| . 
Measurements of pc were made in the two-probe sam- 
ple configuration. Contacts were covering the whole ah 
plane area of the c-axis samples. A four-probe scheme 
was used to measure the resistance down to the contact 
to the sample, i.e. the sum of the actual sample resistance 
Rs and contact resistance Rc was measured. Taking into 
account that Rs » i?c, contact resistance represents a 
minor correction of the order of 1 to 5%. This can be 



resistance represents Re- 

The drawback of the measurement on samples with 
c ^ a is that any inhomogeneity in the contact resis- 
tivity or internal sample connectivity admixes in-plane 
component due to redistribution of the current. To min- 
imize this effect, we performed measurements of pc on at 
least 5 samples of each compositions. In all cases we ob- 
tained qualitatively similar temperature dependences of 
the electrical resistivity, as represented by the ratio of re- 
sistivities at room and low temperatures, pc(0)//9c(300). 
The resistivity value, however, showed a notable scatter 
and at room temperature was typically in the range 1 
to 2 mJlcm. For the sake of comparison we selected the 
samples with the temperature dependence of resistivity 
least similar to that of pa {T) . The value of resistivity for 
these samples at room temperature is shown as a func- 
tion of doping in the top panel of Fig. [T] Typically, these 
samples had the lowest value_of electrical resistivity, as 

This is important since 



described in detail in Ref. 21 



partial exfoliation increases resisitivity values [2l|. As 
a best demonstration of the correctness of our measure- 
ments, thermal conductivity measurements in the normal 
state for samples with a;=0.127, accessed by the applica- 
tion of magnetic field, found Wiedemann-Franz law to 
be obeyed in T — J^O limit (28l |. This quantitative coinsid- 
ence of two independent measurements is very important, 
because cracks can be partially transparent for phonons 
and thus affect thermal and electrical transport in a dif- 
ferent way, leading to gross extrinsic Wiedemann-Franz 
law violation [2^ . The evolutions of the inter- plane resis- 
tivity at room temperature, p(300i4r), and of the residual 
resistivity ratio, pc(Q) / pl^iQQK), with doping are summa- 
rized in Fig.[TJ The resistivity value at room temperature 
for most compositions stays in the range 1 to 1.5 mficm, 
with doping it decreases to approximately 0.5 mficm. For 
several x compositions we were not able to find crystals 
with resistivity values lower than 2 milcm, despite the 
facts that (1) the evolution of the temperature-dependent 
resistivity for these samples followed the general trend, 
(2) close in x compositions show usual resistivity values. 
This limits the accuracy of the absolute pc value deter- 
mination by approximately a factor of two. 

Samples for electrical resistivity measurements with 
current flow along the [100] a-axis in the tetragonal plane 
{pa) were cut into bars of (2 - 3) x (0.1-0.2) x (0.1-0.2) 
mm"^ (a X & X c) . Measurements of pa were made in both 
standard 4-probe and 2-probe configurations and gave 



identical results, see Ref. 26 and 271. Electrical resis- 



directly seen for superconducting samples [2l|, [26|, |27[ at 

temperatures T < Tc, where Rs =0 and the measured 



tivity of the samples at room temperature is shown as a 
function of doping in Fig. [51 Error bar represent statis- 
tisal standard deviation for at least 5 samples of each 
composition. The in-plane resistivity monotonically de- 
creases from 270 pVL.cuY in the parent compound to about 
100 /xJ7.cm in the heavily overdoped composition with 
x=0.48. The magnitude of pa(300) is in good agreement 
with previous report over a narrower doping range [s^]. 
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FIG. 1. Room-temperature inter-plane resistivity, pc{300K), 
of Ba(Fei_a;Coa;)2As2 as a function of doping (top panel). 
Lower panel shows doping dependence of the ratio of re- 
sistivities at low temperatures and at room temperature, 
PciO)/pci300K). 




FIG. 2. Room-temperature in-plane resistivity, pa{SOOK), 
of Ba(Fei-2,Coa;)2As2 as a function of doping (top panel). 
Red stars show resistivity values taken from Ref. 30. Lower 
panel shows doping dependence of the resistivity ratio, 

PaiO)/Pa{300K). 



Residual resistivity ratio shows a rapid increase in the 
range where the Fermi surface topology change (Lifshits 
transition) happens (at x «0.025) [3ll,[32|, reaches a max- 
imum at x=0.05 and then decreases towards minimum 
close to 21=0.1. With further doping the ratio increases, 
the effect which mainly comes from a decrease of resis- 
tivity at room temperature. 

The magnetization measurements were performed on 
cleaved samples to minimize the risk from small amount 
of surface flux. Samples typically had total mass of 10 
to 20 mg. Measurements were performed in a standard 
MPSM SQUID magnetometer in a field of 5 T. Unless 
specially mentioned, magnetization measurements were 
performed in configuration HJ-c. For a composition 
X— 0.325 measurements were also performed with H \\ c. 
They found essentially no anisotropy, similar to our pre- 
vious study 0. 



RESULTS 

In Figs. [3] and U] we present evolution of the 
temperature-dependent resistivity with doping. The 
inter-plane resistivity (top panel. Fig. of the par- 
ent compound decreases sharply below Tsa/j similar to 
the in-plane resistivity (bottom panel. Fig. [3]). In the 
inter-plane resistivity the decrease at Tsm=135 K is pre- 
ceded with resistivity maximum at T* 200 K (shown 
with arrow in Fig. [S]). With doping, the decrease of 
Pc{T) below TsM turns into an increase (as seen for 
samples with a;=0.038 to 0.058), similar to the behav- 
ior of Pair), which shows two anomalies due to split 
structural/magnetic transition [l9| . This change near 
X «0.025 is consistent with the proposed Lifshitz transi- 
tion (Fermi surface topology change) as seen in thermo- 
electric power, Hall effect measurements (sij and ARPES 
[3^ . However, the maximum in Pc{T) at T* remains of 
the same crossover type and does not follow resistivity 
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FIG. 3. Temperature dependence of the inter-plane resistivity, 
pc, normalized by its value at room temperature pc.(?)QQK), for 
samples of Ba(Fei„iCoa;)2 As2 with x < 0.166 (slightly above 
the concentration boundary for the superconducting dome) 
as shown in the figure (top panel). The curves are offset to 
avoid overlapping. Arrows show a position of the resistivity 
ma:ximum, presented as a function of dopant concentration 
in the T — x phase diagram (see Fig. [5] below), cross-arrows 
show position of the resistivity minimum, Tcg, appearing at 
high doping levels. Bottom panel shows doping evolution of 
the temperature-dependent in-plane resistivity, pa, normal- 
ized by room-temperature value pa{300K). Arrows show po- 
sitions of T* and Tcg as determined from pc{T), revealing no 
discernible features in the in-plane resistivity. 



behavior below Tsm (either increase or decrease), sug- 
gesting that it is an independent feature. At doping 
close to optimum, Xopt ~ 0.07, the features due to struc- 
tural/magnetic transition are completely suppressed [in 
both Pa{T) and pc{T)], and the temperature dependence 
of the inter-plane resistivity is dominated by the maxi- 
mum at T* and superconducting transition. 

At the highest doping shown in Fig.[3l a;=0.166, when 
the superconductivity is suppressed, a new feature ap- 
pears in the temperature-dependent interplane resistiv- 
ity: a shallow resistivity minimum appears at Tcg > T* . 
In Fig. m we present the evolution of the resistivity for 
higher Co concentrations, starting from those on the 
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FIG. 4. Temperature dependence of the inter-plane resistivity, 
Pc normalized by its value at room temperature pc(300A'), 
for samples of Ba(Fei_2,Coa;)2As2 with high doping levels 
X >0.127 as shown in the figure (top panel). The curves are 
offset to avoid overlapping. Cross-arrows and stright-arrows 
show positions of the resistivity minimum. Tec, and maxi- 
mum, T*, respectively. Bottom panel shows evolution of the 
temperature-dependent in-plane resistivity, pa, normalized by 
room-temperature value pa{300K). Arrows show positions of 
Tcg and T* as determined from the inter-plane resistivity 
temperature dependence, revealing no discernible features in 
the in-plane resistivity. 



over-doped side of the superconducting dome, x=0.127. 
The top panel shows the inter-plane resistivity, the bot- 
tom panel shows the in-plane resistivity, which shows 
metallic behavior for all compositions. Cross-arrows in 
the top panel show the position of the high-temperature 
cross-over from the metallic to non-metallic temperature 
dependence of the inter-plane resistivity at Tcg- Cross- 
arrows in the bottom panel show the same characteristic 
temperatures with respect to the temperature-dependent 
in-plane resistivity finding no discernible features in the 
Pa{T) curves. 

We summarize the doping evolution of the main fea- 
tures of the temperature-dependent resistivity in the 
phase diagram, Fig. [S] The lines of the superconduct- 
ing, Tc, structural, Ts and magnetic, Tm transitions are 
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FIG. 5. Temperature-doping phase diagram of 

Ba(Fei_a;Coi,)2As2 as determined from inter-plane re- 
sistivity measurements. Inset shows comparison of T* and 
TcG, corresponding to the maxima amd minima in pc{T), 
with TpG as found in NMR study by fiting the temperature- 
dependent Knight shift [SJ, 113] . Arrow in the inset shows 
a minimum estimate for TcG for the border composition 
x=Q.l27. 



discernible in both in-plane and inter-plane resistiv- 
ity. The lines corresponding to maxima, T* , and minima 
TcG of the temperature-dependent inter-plane resistivity 
find no correspondence in the temperature dependence 
of the in-plane transport. 

This phase diagram suggests existence of a critical con- 
centration, at which charge gap vanishes. Interestingly 
enough, at the concentration close to critical, a;cG — 0.30, 
the inter-plane resistivity shows a linear temperature de- 
pendence over a broad temperature range, as seen for a 
sample with a;=0.313 (red curve in the top panel of Fig. 0]) 
for T >20 K. For higher a;, the temperature-dependent 
resistivity develops positive curvature, and can be reason- 
ably described by a sum of T-linear and contibutions, 
Pr.jT) = po + AT + BT^, similar to in-plane transport 
[20l|, as shown in the top panel of Fig. [3] for samples with 
a;=0.343 and a;=0.370. 

In Fig. [6] we compare the inter-plane resistivity with 
earlier evidence of the pseudogap in the electron-doped 
iron arsenides: the temperature dependence of the ^^Co 
NMR Knight shift K and T-normalized NMR relaxation 
rate, as measured in Ref. (2^. We recall that, in 

a simple metal, both K and should be tempera- 
ture independent. In contrast, both Knight shift and 
relaxation rate data in Ba(Fei_a;Coa;)2As2 are strongly 
temperature-dependent. In the parent compound, a;=0. 
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FIG. 6. Comparison of the temperature-dependent in- 
terplane resistivity pc (solid lines, left scale) with the 
temperature-dependent ^''Co NMR Knight shift K{T) (solid 
symbols) and relaxation rate, 1/T\T, (open symbols) from 
Ref. [24| (two right scales) for BaFe2As2 (top left panel) and 
Ba(Fei_2;Coa:)2As2 with x=0.04 (top right panel), x=0.074- 
0.08 (bottom left panel) and x=0. 105-0. 108 (bottom right 
panel). A broad maximum in the temperature dependence 
of the inter-plane resistivity clearly correlates with pseudogap 
features in NMR measurements: a crossover slope change in 
K{T) and the onset of a low-temperature rapid rise in 1/T\T. 



K{T) shows an increase with temperature (seen in all 
compositions), with a mild slope change around 210 K. 
On the other hand, slightly increases on cooling be- 
low 200 K on approaching the temperature of the coupled 
structural-magnetic transition, TgM^ 135 K. These two 
features in K{T) and vs T are close in temperature 
to a shallow maximum in Pc(T) at around 200 K, pre- 
ceding a sharp drop of resistivity at Tsm- 

This correlation between the features in the 
temperature-dependent NMR Knight shift and the 
inter-plane electrical resistivity becomes clearer with 
increasing Co doping. The slope change in the Knight 
shift becomes more pronounced and, for a composition 
with a;=0.105, it shifts down to ~100 K. In both NMR 
measurements and in the inter-plane resistivity the 
features remain of a broad crossover type, with difficult 
to define characteristic temperatures. The resistivity 
maximum is a better defined feature, though it is still 
broad and its location for several samples studied for 
each composition could be slightly affected by the 
admixture of the in-plane resistivity. This admixture 
may affect the Pc{T) for the a;=0.108 sample in Fig.[6l as 
suggested by its slight deviation from the series evolution 
with doping, top panel in Fig. |31 (Small jumps in the 
temperature dependence are extrinsic and are caused by 
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FIG. 7. (Color online) Top panel. Temperature-dependent 
molar magnetic susceptibility, Xmoie{T), measured in mag- 
netic field H ±c= 5 T. The data for low dopings are from 
Ref. [l^. On doping increase the slope of the T-linear increase 
in x{T) (shown with orange line for pure composition x=0 
at T=250 K) decreases and for a;=0.313 the dependence be- 
comes flat at T >~150 K. Bottom panel. Doping evolution 
of the slope of Xmoie{T), dxmoie{T)/dT, at r=250 K. 



sample cracking during thermal cycle.) 

The linearly increasing NMR Knight shift [3j] reflects 
an unusual linear ternperature-dependent static magnetic 
susceptibility, x{T), [35|437|. This anomalous linear x{T) 
dependence was shown to go away at high doping levels, 
being replaced by a Curie- Weiss behavior of susceptibil- 
ity [35| . The magnitude of the Knight shift variation also 
diminishes for over-doped compositions, and it was sug- 
gested that the pseudo-gap feature disappears at critical 



33| 



concentration for superconductivity xsc ~0.2 

We performed magnetization measurements on Co 
doped compositions, as shown in Fig. [T] The behavior 
at low Co-dopings < x < 0.166 was studied system- 
atically by Ni et al. [3| for both orientations of mag- 
netic field parallel and perpendicular to the tetragonal 
c-axis and found small anisotropy. Our new measure- 
ments for X >0.166 concentrate on doping evolution of 



susceptibility in configuration HJ^c. We performed mea- 
surements with H \\ c only for sample a;=0.343 and found 
no anisotropy. As can be seen from Fig. [3 the slope of 
the T-linear portion in x(^) gradually diminishes with 
X. The Xmoie{T) curve for a;=0.290 shows very small 
but still clearly discernible increase with T, though with 
the slope notably smaller than the slope for a;=0.270. For 
a;=0.313 the increase of magnetic susceptibility with tem- 
perature is completely gone. Instead, x(^) becomes tem- 
perature independent above 150 K. The Curie- Weiss in- 
crease of x{T) on cooling at low temperatures for sample 
with x=0.313 is most likely extrinsic, it is not observed 
for samples with lower and higher doping, x=0.290 and 
a;=0.343. On the other hand, the XmoieiT) dependence 
does not reveal any increase at high T for both a;=0.313 
and two different samples of a;=0.343. 

To quantify this evolution, in the bottom panel of 
Fig. [7] we show a slope of Xmoie{T) at r=250 K as a 
function of x. The dependence shows a dramatic change 
between x=0.290 and a;=0.313, the same concentration 
as xcG determined from inter-plane resistivity. We note 
that since the crystals used in the inter-plane resistiv- 
ity and magnetization measurements are from the same 
batches, there is minimal uncertainty in the concentra- 
tion comparison. Both inter-plane resistivity and mag- 
netic susceptibility show pronounced changes of behavior 
between a:=0.290 and a;=0.313. For x >0.313, the inter- 
plane resistivity increases monotonically and superlin- 
early with temperature as expected for a metal. The flat 
temperature-independent magnetic susceptibility is ex- 
pected for Pauli susceptibility of a metal as well. The dis- 
appearence of the linearly rising x(T) for heavily doped 
compositions is similar to the early report [s^, though 
with notable difference in the concentration boundaries 
(x=0.125 in our notations vs 0.30). This discrepancy 
may be a result of poor composition control in early 
crystals. We would like to point out that at similar con- 
centration ~0.4 (there is no systematic doping-evolution 
study in the range) Hall constant becomes temperature- 
independent, again in line with expectation of the behav- 
ior of a usual metal Isa . 



DISCUSSION 
Doping evolution of the anisotropy 

The electronic structure of Co-doped BaFe2As2 is now 
well established to be three-diemnsional by various tech- 
niques 3^43|. However, evolution of the anisotropy with 
doping was never studied in a systematic way. From 



Fig. [2] we can see that doping in the range from x—Q 
to a;=0.48 leads to an approximately 3 times decrease of 
the in-plane resistivity at room temperature, agre eing , 
within error bars, with the previous measurements [30|. 
For pc, Fig.[TJ the variation is approximately of the same 
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magnitude, keeping in mind an uncertainty of the factor 
of 2 for pc values. This suggests that the evolution of 
the anisotropy in a broad doping range is very gradual, 
and puts an upper bound of about a factor of two for the 
anisotropy change. Comparison of these numbers with 
the existing band structure calculations [2l|, |22|, 144447 
should be taken with a grain of salt, since variation of 
the position of the As atom in the lattice from the one 
obtained in experiment to the one calculated from total 
lattice energy minimum 2l|, |22| brings the effect which by 
far exceeds the total anisotropy variation. A general de- 
crease of both in-plane and inter-plane resistivities with 
doping is suggestive that charge is actually donated into 
the system, which does not go in line with suggestions 
that all doping can be treated as additional scattering 
i3|- 

The anisotropy at low temperatures, important for the 
anisotropy of the upper critical field in the superconduct- 
ing state, is heavily affected by the structural/magnetic 
transition and by the pseudogap. We will discuss these 
effects in the next sections. 



Structural/magnetic ordering and inter-plane 
resistivity 

Magnetic ordering below Tsm , presumably of spin den- 
sity wave (SDW) nature [il], reconstructs the Fermi sur- 
face, opening a superzone gap on electron and hole pock- 
ets. This is seen as an increase of the in-plane resistivity 
in BaCol22 with x=0.025 to 0.058. The parts of the 
Fermi surface which are not affected by the superzone 
(SDW) gap, enjoy a notably reduced inelastic scattering 
in the magnetically ordered phase (g^HO,!!^. In the par- 
ent compound, in which elastic scattering is small, this 
decrease of scattering overcomes the loss of the carrier 
density so that the total conductivity increases below 
Tsm- Since the inter-plane transport is determined by 
the warped parts of the Fermi surface [2l| , least affected 



by the SDW superzone gap, the inter-plane resistivity 
should be affected much less by the SDW gap opening 
than Pa- This is indeed seen, in BaCol22, most clearly 
for sample with a;=0.012. Here, the in-plane resistivity 
shows an intermediate behavior between pure and heav- 
ier doped compositions: whereas Pa{T) increases imme- 
diately below Tsm and then shows a shallow decrease to 
a much higher residual value than in pure samples, the 
inter-plane resistivity does not manifest a local maximum 
below Tsm and the resistivity decrease is almost as steep 
as in pure crystals. The features in the temperature- 
dependent resistivity upon crossing structural and mag- 
netic transitions [19|, l50| can be similarly resolved in in- 
plane and inter-plane transport, though structural tran- 
sition is always less pronounsed in Pc{T). 



Maximum in the temperature-dependent inter-plane 
resistivity at T* 



The decrease of the inter-plane resistivity below T* 
shows a clear correlation with the NMR Knight shift, 
therefore we need to look for a common origin. An 
important observation is that in BaCol22 with a;=0.08. 
Fig. ini both the Knight shift and the inter-plane resistiv- 
ity at low temperatures, T < T*, follow expectations of 
a metal with temperature-independent density of states: 
the resistivity shows metallic decrease on cooling, and 
the Knight shift is temperature-independent. Simultane- 
ously, the increase indicates slowing down of mag- 
netic fluctuations. This suggests that magnetic correla- 
tions may play an important role in the anomalies in all 
three measurements. Same trend hold for sample with 
x=0.105, however, the features in NMR measurements 
fade away with overdoping. 



In the NMR study of Ref. 17, the temperature- 
dependent Knight shift was fit using a two-component 
model, with K = A + B * exp{-TpG/T). At tem- 
peratures T <C TpG this crosses-over to a metallic be- 
havior with constant Knight shift detemined by the A 
term unaffected by the pseudo-gap. At T > Tpg both 
terms become temperature-independent and we can ex- 
pect restoration of the metallic behavior with K — A+B. 
Fitting the temperature-dependent Knight shift, the au- 
thors determined rpG=560 K ± 150 K for optimally 
doped BaCol22 samples, a;=0.08 (with A=0.715% and 
B=0.244%), Ref. M, and rpG=450 K for a;=0.26 (with 
A=0.20% and B=0.23%). Assuming that the A and B 
coefficients represent partial DOS contributions of the 
ungapped and gapped parts of the Fermi surface, respec- 
tively, we would expect that at temperatures of the order 
of Tpg/ 3 or so, which would give us a temperature in 
100 to 200 K range for optimal doping, we restore metal- 
lic resistivity temperature dependence, while resistivity 
decrease with temperature would be observed at higher 
temperatures due to carrier activation. This is consistent 
within general trend in evolution of pdT), but not Pa{T)- 
This would suggest that the pseudogap affects predomi- 
nantly the most warped parts of the Fermi surface. 



We need to notice though, that it is difficult to explain 
resistivity decrease above T* solely by the existence of 
a spin gap, as probed by NMR [3J]. Activation of spin 
fluctuations in the metallic phase can only increase scat- 
tering of charge carriers, which is seen in in-plane trans- 
port. Decrease of the inter-plane resistivity, despite being 
very small, would require rather an increase of the carrier 
density by excitations over a charge gap. 
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Minimum in the temperature-dependent 
inter-plane resistivity and charge gap 



The importance of charge gap formation for non- 
metalhc temperature dependence of pc above T* is most 
clearly seen from the temperature-dependent inter-plane 
resistivity for BaCol22 composition with 0.166 < x < 
0.270. Here the resistivity shows metallic temperature 
dependence at high temperatures and then crosses over to 
a non-metallic increase below a temperature of a charge 
gap formation Tqq. The monotonic evolution of the 
curves suggests that for lower dopings, x <0.166, the 
TcG goes above the experimentally accessible range. If 
this is true, the end of the temperature range of mono- 
tonic resistivity decrease on heating gives us an estimate 
for a minimum value of charge gap for compositions with 
X < 0.127 as above 400 K. In inset in Fig. [5l we compare 
experimentally determined TcG from the inter-plane re- 
sistivity measurements with Tpo determined from fitting 
K(T) curves. Both measurements have very big error 
bars, and yet they do not match well. This may suggest 
another possibility, that a metallic temperature depen- 
dence of Pc at high temperatures is only confined in some 
range of dopings. 

The scenario with the existence of a semi-metallic 
charge gap was invoked for an explanation of the T- 
linear magnetic susceptibility, with simultaneous str ong 
teniperature-dependcnt Hall and Seebeck coefhcients [30j, 
3ll . l5ll | . In this model thermal activation of carriers over 
a narrow gap results in a carrier density increase with 
temperature. This would naturally lead to a decrease of 
the inter-plane resistivity with temperature. We should 
notice though that the magnitudes of the effects, nec- 
essary to explain temperature-dependent magnetization, 
by far exceed the magniture of the pc variation observed 
in our experiments. This is also true with respect to dop- 
ing evolution of the characteristic temperatures, Tpc and 
TcG, Fig- El which do not connect gradually. Simultane- 
ously the linear rise in magnetization with temperature 
does not coincide with resistivity maximum in our study, 
especially for pure BaFe2As2. 



Despite this clear discrepancy between the two sug- 
gested explanations for the temperature-dependent mag- 
netization and transport and our data, the effects in 
the inter-plane resistivity and in the magnetization are 
clearly correlated. In addtion, determination of the char- 
acteristic pseudogap temperature Tpc in NMR measure- 
ments is heavily model dependent, whereas the minimum 
in the temperature dependence of the inter-plane resistiv- 
ity, despite being boad, is rather well defined for x >0.166 
and shows a systematic evolution. This may imply that 
we need to invoke different mechanisms for the explana- 
tion of the pseudogap features. 



Critical concentration 

We now turn to the evolution of Pc(T) and magnetic 
susceptibility in the vicinity of xcg- The most remark- 
able observation here is that at a;=0.313 the resistivity 
is faily linear over a broad temperature range from ap- 
proximately 400 K down to 20 K and saturates at lower 
temperatures. For a;=0.290 the dependence is also close 
to T-linear with a shallow slope change at '^150 K. For 
doping with x > 0.313 the Pc{T) becomes superlinear, 
similar to pa(T), and its inelastic part can be reasonably 
decribed as a sum of T-linear and terms, as shown 
for the curves a;=0.343 and x=0.370. In general, evolu- 
tion of Pc(T) with doping is reminescent of the one found 
in systems on the verge of magnetic order and assigned 
to the existence of the magnetic quantum critical point. 
This observation suggests that the pseudogap is magnetic 
in origin, and is accompanied by the charge gap, rather 
than the charge gap itself being responsible for anoma- 
lous electronic properties. 

On note, none of the anomalies in the magnetic prop- 
erties is clearly reflected in the in-plane transport. This 
unusual single-axis effect of the pseudogap on the resis- 
tivity suggests that the magnetic action is concentrated 
on a small fraction of the Fermi surface, and importantly, 
on the most warped part contributing mainly to the inter- 
plane transport. 



Origin of the pseudogap 

The existence of two additional crossover lines in the 
phase diagram of Ba(Fei_a;Coa;)2As2, as revealed by 
the inter-plane resistivity, raises an interesting question 
about their origins. Strong anisotropy of the pseudogap 
makes a scenario, in which the gap is due to supercon- 
ducting pairing of charge carriers, however without su- 
perconducting condensate formation, very unlikely. We 
therefore should consider the possibility that the pseu- 
dogap is arising due to either short range, or short-lived, 
magnetic correlations, or represents a partial gap in the 
electronic structure. 

The magnetic structure of parent BaFe2As2 is char- 
acterized by a stripe-type antiferromagnetic ordering, in 
which antiferromagnetic spin arrangement is typical for 
directions both in the plane and between the planes, 
introducing three-dimensional magnetic Brillouin zone, 
poorly matching the Fermi surface. It is difficult to ex- 
pect pronounced anisotropy for this case. On the other 
hand, if correlations seen by the inter-plane transport 
were the same as those of the ordered phase, it would be 
difhcult to explain a pre-transition decrease of resistivity 
below T* in Co-doped samples with a:=0.037 to 0.058, 
with successive increase of resistivity below Tsm- This 
may suggest that uniaxial anisotropy of the pseudo-gap 
comes from magnetic fluctuations with a different char- 
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acteristic wavevector. A situation like this, when fluc- 
tuations and ordering wavevector are not the same was 
found in some intermetallic and heavy fermion systems 
(5^ . Indirect evidence for such a possibihty comes from 
the fact that in a closely related EuFe2As2, antiferro- 
magnetic ordering of Eu moments happens between the 
planes, while the Fe layer moments remain parallel in the 
planes 5^, I^l • Since this ordering is mediated by RKK Y 
interactions via conduction electrons, it suggests that the 
generalized spin susceptibility may have maxima, which 
correspond to the existence of the inter-plane nesting. 

In closely related BaMn2As2, magnetism is of local mo- 
ment type, and the magnetic order is of usual AF G-type 
[ssf . In EuRh2As2 , commensurate and incommensurate 
spiral-like structures with propagation along the c-axis 
are found. 56| Although these compounds are differing 
in band structure and electron count, these observations 
of different types of ordering may be suggesting that var- 
ious magnetic structures are not very different in energy. 

In the lack of any evidence for the existence of such 
correlations in Co-doped BaFe2As2, we just speculate 
what consequences uniaxial character of the pseudogap 
may have. This type of a pseudogap is impossible in 
two dimensional cuprates, it is a direct consequence of 
the difference in the dimensionality of the electronic and 
magnetic systems in the cuprates and in iron arsenides. 
If the link between the symmetry of the pseudogap and 
of the superconducting order parameter, as found in the 
cuprates 0,[53|, is preserved in the iron arsenides, c-axis 
pseudogap would correspond to a superconducting gap 
having maxima/minima at the poles. This scenario was 
invoked theoretically for explanation of unusual behavior 
in the superconducting gap [s^. In experiment, varia- 
tion of the superconducting gap with polar angle is found 
in inelastic neutron scattering in Ni doped compound 
at optimal doping [soj, with gap magnitude decreasing 
towards the poles, and in penetration depth study of 
BaNil22 [g^. It is directly revealed in the inter-plane 
heat transport study 28], as opposed to the in-plane 
study jell. 

Finally we would like to point to a certain similarity 
in the critical behaviour of the inter-plane resistivity in 
BaCol22 and in CeCoIns. In CeCoIn5, a true critical 
behavior at a field-tuned QCP 62, 6^ with T- linear re- 
sistivity and violation of the Wiedemann-Franz law is 
observed for transport along the tetragonal c-axis [63 |. 
Transport in the plane perpendicular to c-axis, despite 
showing unusual power law behavior, obeys the WF law 

M. 



inter-plane resistivity, pc{T), reveals anomalous fea- 
tures clearly correlating with features in the tempera- 
ture dependence of the the NMR Knight shift and spin- 
relaxation rate, assigned to the formation of the pseudo- 
gap. Evolution of Pc{T) with doping reveals two char- 
acteristic energy scales, of the the resistivity maximum 
(seen for compositions < x <^ 0.2) and resistivity min- 
imum at a tempearture Tqq, seen for 0.166 < x < Xc, 
Xc «0.3. The temperature-dependent pc is close to lin- 
ear close to a;cG and super- linear for x > xqg- None 
of these features are evident in the in-plane resistivity 
Pa{T). For doping levels x < xqg, x{T) shows a known, 
anomalous, T-linear dependence, which disappears for 
X > xcG- These features are consistent with the exis- 
tence of a uniaxial charge gap, accompanying formation 
of the magnetic pseudogap, and its critical suppression 
with doping. This evolution suggests existence of critical 
point for pseudogap order. The superconducting dome is 
confined inside the pseudo-gap dome. 
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CONCLUSION 

Contrary to the in-plane electrical resistivity, which 
away from the domain of structural/magnetic order- 
ing shows monotonic metallic temperature dependence. 
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